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THEORY  OF  ELECTRON  CYCLOTRON  RESONANCE  HEATING 

OF  TOKAMAK  PLASMAS 


I.  INTRODUCTION 

Radio  frequency  plasma  heating  is  of  great  potential  importance  in  controlled  fusion 
research.  Heating  at  the  electron  cyclotron  frequency  is  one  form  of  plasma  heating  which 
offers  many  potential  advantages.  Until  recently,  however,  electron  cyclotron  resonance  heat¬ 
ing  has  received  relatively  little  attention,  primarily  because  efficient  powerful  microwave 
sources  in  the  millimeter  wavelength  range  have  not  been  available.  This  situation  is  now 
changing  due  to  the  development  of  the  "gyrotron"  device  both  in  this  country  and  in  the  Soviet 
Union  The  gyrotron  produces  microwave  power  by  the  interaction  of  a  high  energy  anisotro¬ 
pic  electron  beam  with  electromagnetic  waves  via  the  cyclotron  resonance.  In  contrast  with 
conventional  tubes,  no  cavity  resonance  is  involved,  and  thus  much  shorter  wavelengths  can,  in 
principle,  be  obtained.  Initial  tokamak  plasma  heating  experiments  at  the  electron  cyclotron 

frequency  have  now  been  carried  out,  and  more  are  planned. 

>■- 

Motivated  by  these  developments  we  (and  others^)  have  attempted  a  theoretical  study  of 
tokamak  plasma  heating  by  externally  launched  electromagnetic  waves  whose  frequency 
matches  the  electon  cyclotron  frequency  in  the  interior  of  the  device.  Section  II  reviews  the 
basic  limitations  imposed  by  wave  accessibility  to  the  center  of  the  tokamak.  It  is  shown  that  a 
density  limit  for  accessibility  exists.  In  Sec.  ill  we  utilize  the  ray  equations  to  obtain  the  wave 
propagation  and  energy  deposition  in  a  self-consistent  tokamak  equilibrium  (this  represents  the 
first  such  study  utilizing  a  real  toroidal  plasma  equilibrium  and  including  important  relativistic 
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effects  on  the  wave  absorption).  We  find  several  favorable  features  of  this  heating  method. 
Sections  IV  and  V  investigate  two  possible  mechanisms  which  might  substantially  alter  the 
favorable  picture  obtained  in  Sec.  Ill:  scattering  of  the  incident  heating  wave  by  low  frequency 
density  fluctuations  existing  in  the  tokamak  (e.g.  drift  waves),  and  the  excitation  of  parametric 
instabilities  by  the  heating  wave.  Our  investigation  suggests  that  these  two  mechanisms  should 
not,  in  fact,  pose  a  serious  problem  for  electron  cyclotron  heating  of  tokamak  plasmas.  In  par¬ 
ticular,  the  threshold  incident  power  density  for  the  excitation  of  parametric  instabilities  should 
far  exceed  the  actual  incident  power  density  to  be  expected  in  typical  situations.  Section  VI 
summarizes  results  and  presents  conclusions. 

II.  REVIEW  OF  ACCESSIBILITY 

For  simplicity,  in  this  section  we  consider  a  straight  plasma  cylinder  with  density 
N(r)(dN/dr<0  and  N(a)=0),  and  magnetic  field  B  —  B(r,0)zo.  To  simulate  the  toroidal  va¬ 
cuum  field  dependence  on  the  m^jor  radius  we  take  B—B„(\-x/R„)~ 1  where  x— rcos0  and 
R„  is  the  major  radius  of  the  torus.  Furthermore,  we  .assume  incident  waves  launched  at  0-0 
or  0-ir  with  no  0  component  of  the  incident  wavenumber  (fc,-0).  Since  B  has  no  0  com¬ 
ponent,  k»  remains  zero  as  the  wave  propagates  inward.  (Note  that  these  simplifications  are  not 
made  in  our  ray  tracing  calculations  to  be  presented  in  Sec.  III.) 

A.  Normal  Incidence 

a.  Ordinary  mode:  The  dispersion  relation  is 

ij2  -  \-X 

where  tj2  -  Ar2c2/«2,  X  -  «2/w2,  un  is  the  plasma  frequency,  k  is  the  wavevector,  c  is  the 
speed  of  light  and  o»  is  the  wave  frequency.  To  avoid  reflection  we  require  X„  <  1  (Fig.  I), 
where  X„  is  X at  r-0.  For  ai-oi,®  the  condition  is  Wp,,  <  «(®,  where  and  «,°r  are  the  plasma 
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frequency  and  the  electron  cyclotron  frequency  («„)  evaluated  at  /--0.  [For  first  harmonic 
heating,  g>-2&>°,,  the  condition  becomes  to^  <  2 and  heating  of  higher  density  plasmas  be¬ 
comes  possible.  First  harmonic  heating  will,  however,  not  be  considered  in  this  paper.] 

b.  Extraordinary  mode:  The  dispersion  relation  is 

v2-l-[X(l-X)/(l-X-Y2)] 

where  Y  -  <oce/to.  The  zeros  of  ij2  (turning  points)  are  located  at  X  -  1  ±  Y.  A  pole  of  tj2 
(resonance)  is  located  at  X  —  l-Y2  and  is  the  condition  for  upper  hybrid  resonance.  For 
to  —  a>ce,  Y-(l-x/R0)-'  and  i)2  as  a  function  of  x  is  plotted  in  Fig.  2(a)  for  X0  <T  2.  We  see 
that  an  extraordinary  wave  launched  from  the  low  magnetic  field  side  (x— a)  always  en¬ 
counters  a  turning  point  (followed  by  the  upper  hybrid  resonance)  near  the  plasma  edge.  Thus 
the  extraordinary  mode  should  be  launched  from  the  high  magnetic  field  side  where  it  will  pass 
through  the  cyclotron  resonant  absorption  region  near  x-0.  If  X0>2,  the  wave  will  be 
reflected  for  both  low  and  high  field  side  injections  (Fig.  2(b)). 


where  q  =  ckjto,  S-sin&-k„c/<o,  0  is  the  angle  of  incidence  (we  have  looked  at  normal  in¬ 
cidence  5-0  in  Sec.  A),  and  ^  and  *n  denote  components  of  k  perpendicular  and  parallel  to  B. 
The  plus  sign  in  front  of  the  square  root  is  used  for  the  ordinary  mode,  and  the  minus  sign  for 
the  extraordinary  mode  (this  can  be  seen  by  letting  S2-*0  in  (1).) 
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>  1  + 


1  — S2 
4 S2  ’ 


and  Xj  <  1,  then  the  wave  will  be  mode  converted  to  ordinary  wave  before  reaching  the  upper 
hybrid  resonance.  The  ordinary  wave  is  reflected  at  X— X\,  moves  back  to  the  mode  conversion 
point,  and  is  converted  back  to  extraordinary  wave  which  then  propagates  back  towards  the  high 
held  side  (  Fig.  3  ). 


In  electron  cyclotron  resonance  heating  experiments,  we  want  to  avoid  all  reflection,  reso¬ 
nance,  and  mode  conversion  points  before  reaching  the  center  of  the  torus  where  maximum 
absorption  occurs.  By  examination  of  the  expressions  for  Xu  X2.Xy  Xr  and  Xc  we  find  that  this 

.  ,  .  '  *  ‘  •  .  t 

can  be  achieved  for  K„- 1  if 


X0 


2  (1-S2),  for  S2  ^  1/3. 

1+  (1~522)2  ,  for  S2S  1/3, 
4«S 


(8) 


for  the  extraordinary  wave  mode  launched  from  the  high  magnetic  field  side,  and  if 

jl,  for  S2^  1/2, 

X°  <  |2(1-S2),  for  S2>  1/2,  (9) 

for  the  ordinary  wave  mode  (launched  from  either  side).  These  accessibility  conditions  are 
displayed  graphically  in  Fig.  4.  [Note  that  the  extraordinary  wave  with  Y0  —  1  is  never  accessi¬ 
ble  when  launched  from  the  low  magnetic  field  side  since  it  will  always  encounter  the  reflection 
X“Xj  close  to  the  plasma  edge  and  the  resonance  X—Xr  soon  after  (as  in  Fig.  2).) 


In  order  to  discuss  the  accessibility  problem  for  a  reactor  plasma  it  is  instructive  to 
express  ot}/<o2e  in  terms  of  /9,  the  ratio  of  plasma  to  magnetic  field  pressure, 

w2/«2,  -  1.25  (/3/10%)  (10  keV/r,). 
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where  we  have  assumed  equal  electron  and  ion  temperatures,  Tr  -  Tr  It  is  currently  thought 
that  an  upper  limit  on  /3  is  set  by  the  stability  of  the  plasma  to  magnetohydrodynamic  balloon¬ 
ing  modes.  For  example,  if  the  central  /3  is  limited  to  below  16%  at  Tf  -  10  keV,  then  the 
extraordinary  mode  is  necessarily  accessible,  while  for  an  8%  limit  both  the  ordinary  and  the 
extraordinary  modes  are  accessible. 

III.  LINEAR  THEORY  OF  WAVE  PROPAGATION 
AND  ENERGY  DEPOSITION 

In  this  section  we  utilize  the  geometrical  optics  approximation  to  study  the  propagation 
and  energy  deposition  of  ordinary  and  extraordinary  mode  waves  in  a  tokamak  plasma.  Our 
formulation  will  be  valid  if  the  plasma  inhomogenity  scale  length  and  the  damping  length  of  the 
wave  greatly  exceed  the  wavelength  of  the  wave.  These  conditions  are  well-satisfied  in  all  the 
cases  studied.  In  Sec.  Ill  (a)  we  discuss  the  formulation  of  the  ray  equations  in  toroidal  coordi¬ 
nates,  the  local  dispersion  and  damping  of  the  waves,  and  the  self-consistent  tokamak  equilibri¬ 
um  used.  Section  III  (b)  presents  and  discusses  the  results  of  computer  calculations  based  on 
the  equations  developed  in  Sec.  Ill  (a)  and  Appendix  A. 

A.  Formulation 


We  consider  waves  satisfying  a  local  dispersion  relation  D(«  -  i  y,  r,  k)  -  0,  where  <v 
and  y  are  the  real  frequency  and  damping  rate  of  the  wave  and  k  is  the  wavenumber.  Dean  be 
separated  into  two  parts  D  -  D,  +  /  D,  where  D,  and  D,  are  real  for  y  -  0.  Then  for  y  «  m 
we  can  approximate  w  as  the  solution  of  D,(o>,  r,  k)  -  0  and  y  as 
y  s  T(r,  k)  -  D,(«,  r.  k)  (dD,/d«)-1.  In  the  geometrical  optics  limit  the  orbit  trajectory  and 
damping  of  the  energy  in  a  wave  packet  are  then  given  by  the  ray  equations.’ 


dk  _ 

dD, 

dt 

dD, 

dr 

Or 

'dr 

dk  ’ 

dt_  . 

d  Dr 

dw 

-  -2T w, 

dr 

do* 

dt 

6 

(10a, b) 
(10c,d) 
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where  w  is  the  energy  in  the  wave  packet.  The  quantity  t  is  a  measure  of  the  distance  along 
the  ray  trajectory.  Equations  (10b)  and  (10c)  yield  the  usual  group  velocity: 
dt/dt  —  —(bDr/bk)/(bDr/bu>)  =  9a>/dk.  Note  that  in  Cartesian  coordinates  Eqs.  (10a, b)  are 
identical  to  Hamilton’s  equations  of  motion  for  a  particle  with  momentum  k  and  Hamiltonian 
D„  dkjdx  —  -  bDJbXj  and  dxjdi  «  bDJbkj.  We  now  introduce  the  usual  toroidal  coordi¬ 
nates  (r,  9,  <t>),  where  r  is  the  minor  radial  position,  9  is  the  poloidal  angle,  and  </>  is  the 
toroidal  angle.  The  distance  from  the  major  axis  of  the  torus  to  a  point  (/-,  9,  <t>)  is 
/?-/?„  +  rcosfl.  It  is  convenient  to  preserve  the  Hamiltonian  form  of  (10a,  b)  which  exists 
in  Cartesian  geometry  by  utilizing  momenta  variables  which  are  cannonically  conjugate  to 
(r,  9,  <f>).  Therefore  we  introduce  kr,  m  —  /•  A»  and  n  —  R  A^  where  (Ar,  k„,  k ^)  are  the 
( r ,  9,  <t> )  components  of  k.  In  terms  of  (r,  9,  <t>)  and  their  connonically  conjugate  "momenta" 
( kr ,  m.  n)  the  ray  equations  (10a, b)  become6  7 

dr  _  d£  _  0^  d±  =  bDr 

dr  dkr  ’  dr  dm  dr  bn 

dkr  ^  bD,  dm  9Dr  drt  _  „  .  . 

dr  br  dr  b9  '  dr 

and  we  have  assumed  toroidal  symmetry  so  that  bDr/b<t>  =  0. 

Since  the  wave  phase  velocity  is  much  larger  than  the  electron  thermal  speed  we  may  use 
the  cold  plasma  dispersion  relation  for  D„  Eq.  (1).  In  (1)  we  now  use  for  oj,,  and  a>,,.  their  lo¬ 
cal  (r-dependent)  values,  set  S  -  k„c/io,  AM  -  B  •  k/B  -  ( B,k,  +  B„m/r  +  B6  n/R)B~ 1  with 
B(r)  the  equilibrium  tokamak  magnetic  field,  and  k l  *■  A2  -  A,?  — 

A,2  +  ( m/r )2  +  («//?)2  -  A2.  The  damping  rate  T  is  due  to  electron  cyclotron  resonance  of  the 
wave  with  particles  satisfying  the  condition 

»  "  Mu  +  w«-/y» 

where  y„  is  he  relativist  mass  factor,  y„  -  [1  -  (v/c)2]_1/2,  v  is  the  particle  velocity  and  vM  — 
v  ■  B /B.  The  calculation  of  T  is  straight-forward  and  is  given  in  Appendix  A. 
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To  complete  the  formulation  we  need  to  specify  the  equilbrium  magnetic  field,  B,  plasma 
density,  N ,  and  electron  temperature,  Te  as  functions  of  the  coordinates  r  and  0.  (We  assume 
Te  —  T,  subsequently.)  These  are  then  substituted  in  the  expressions  for  D,  and  T  to  complete¬ 
ly  specify  their  spatial  dependence.  We  assume  that  the  toroidal  component  of  B  is  a  vacuum 
field,  B*  ~  R~\ 

B+  -  B0[  1  +  cos0(r/Bo)]-1. 

The  other  equilibrium  quantities  are  obtained  by  considering  a  particular  straight  cylinder  equili¬ 
brium  and  then  making  use  of  Shafranov’s  aspect  ratio  expansion*  to  obtain  toroidal  corrections 
to  this  equilibrium.  Essentially,  this  gives  a  self-consistent  toroidal  equilibrium  which  is  correct 
up  to  (and  including)  terms  of  order  ( r/Rn ).  For  the  straight  cylinder  we  choose 

N(r)  —  N{o)  ll  ~{rla)1) 

B„(r)  -  C,(r/R„)  B„[\  + 

where  r  —  a  is  the  plasma  edge,  and  C\  and  C 2  are  constants.  The  details  of  the  toroidal  bend¬ 
ing  of  the  straight  cylinder  with  these  density  and  B„  profiles  are  given  in  Ref.  6  and  will  not  be 
repeated  here. 

B.  Results  and  Dicusslons 

We  have  calculated  the  absorption  properties  of  three  different  types  of  tokamaks:  (A) 
"small  size",  an  example  is  the  Versator  II  (a  -  11  cm,  R„  -  40  cm)  at  the  Massachusetts  In¬ 
stitute  of  Technology;  (B)  "present-day  size",  an  example  is  the  Princeton  Large  Torus  (a  -  40 
cm,  R„  -  132  cm);  and  (C)  "reactor  size",  an  example  is  the  UWMAK  111  (a  —  270  cm,  R„  - 
810  cm),  a  conceptual  reactor  design  done  at  the  University  of  Wisconsin.9 

In  our  calculations,  the  ordinary  mode  is  always  launched  from  the  weak  magnetic  field 
side,  while  the  extraodinary  mode,  is  launched  from  the  high  field  side  (to  avoid  the  cutofT). 
We  also  usually  eschew  launching  the  extraordinary  mode  at  0„  -  180*  (where  9n  is  the 
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poloidal  angle  at  which  the  wave  is  launched,  cf.  Fig.  5)  because  the  inside  of  the  torus  is  usu¬ 
ally  occupied  by  transformers  and  other  hardware,  and,  consequently,  there  is  not  much  physi¬ 
cal  space  available  for  wave  launching  structures.  Instead,  we  shall  usually  take  0o  ~  130°. 

For  the  extraordinary  mode,  the  ray  equations  are  solved  until  the  ray  reaches  the  upper 
hybrid  resonance.  At  this  point,  the  cold  plasma  approximation  for  Dr  breaks  down,  and  the 
calculation  is  stopped.  As  shown  by  Stix,10  the  extraordinary  wave  mode  converts  to  a  Bern¬ 
stein  wave  at  the  upper  hybrid  resonance.  The  Bernstein  wave  will  then  propagate  back  into 
the  region  of  increasing  magnetic  field.  We  have  examined  the  propagation  and  absorption  of 
Bernstein  waves  originating  from  the  upper  hybird  resonance,  and  we  have  found  that  the 
waves  propagate  back  toward  the  center  of  the  torus  and  are  typically  fully  absorbed  shortly  be¬ 
fore  they  reach  the  cyclotron  resonance. 

The  calculations  which  we  present  here  are  of  most  use  when  the  launched  ordinary  or 
extraordinary  wave  is  absorbed  essentially  completely  in  one  pass  through  the  center.  As  we 
shall  see,  this  is  typically  the  case  for  the  larger,  hotter  tokamaks  (e  g.  Princeton  Large  Torus  to 
reactor  size),  but  not  for  smaller  devices  (e.g.  Versator  ID.  For  the  smaller  devices,  the  full 
picture  of  extraordinary  mode  absorption  involves  the  Bernstein  wave  mode  conversion  process 
discussed  previously;  while  for  the  ordinary  mode,  reflections  from  the  tokamak  walls,  and  the 
limiter,  etc.,  must  be  considered. 

Type  A  Tokamaks 

In  Table  I,  we  show  the  absorption  of  ordinary  and  extraordinary  modes  versus  angle, 
density,  temperature,  and  magnetic  field  for  type  A  tokamaks  (i.e.  small  size  tokamaks).  The 
definitions  of  symbols  appearing  in  this  table  are  as  follows.  0„  denotes  the  poloidal  angle  loca¬ 
tion  at  which  the  wave  is  launched.  (  is  the  angle  of  incidence  in  the  poloidal  plane;  i.e.,  the 
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angle  between  the  projection  of  the  launched  wavenumber  on  the  poloidal  plane  and  a  line  join¬ 
ing  the  launch  point  with  the  point  /•— 0.  The  definitions  of  the  angles  9„  and  (  are  shown  in 
Fig.  5.  The  angle  9  specifies  the  ray  direction  at  the  launch  point  with  respect  to  the  toroidal 
direction, 

sin  9  -  (n/f?)/(w/c) 

where  R  is  evaluated  at  the  launch  point.  rm„  is  the  radius  of  the  magnetic  surface  where  the 
energy  deposition  per  unit  volume  averaged  over  a  magnetic  surface  is  maximum.  Also,  the 
symbols  X  and  O  in  the  table  denote  extraordinary  and  ordinary  modes.  We  see  from  Table  I 
that,  for  this  type  of  tokamak,  the  absorption,  in  general,  is  small  (due  to  the  relatively  low 
temperatures  involved).  While  the  density  dependence  for  absorption  is  weak  for  the  ordinary 
mode,  it  is  strong  for  extraordinary  mode.  In  fact,  the  extraordinary  wave  attenuation  may  be 
shown  to  be  inversely  proportional  to  density.  The  absorption  for  both  modes  does  not  vary 
much  over  a  wide  range  of  9„,  except  near  9„  —  n/2.  If  the  angle  9„  approaches  ir/2  (i.e.,  the 
wave  is  launched  close  to  the  resonance),  the  energy  is  absorbed  near  the  edge. 

The  angular  dependence  on  9  (which  is  a  measure  k\\)  for  the  ordinary  mode  is  also  weak 
near  9  —  0.  The  ordinary  wave  absorption  decreases  slowly  as  9  increases.  In  contrast,  the 
extraordinary  mode  absorption  depends  strongly  on  9,  and  increases  as  9  increases,  until  the 
mode  encounters  the  turning  point  X2.  We  note  that  (from  the  small  values  of  rmJa  in  the 
table)  for  both  modes  most  of  the  energy  absorbed  is  deposited  near  the  center. 

The  transmission  profiles  for  the  ordinary  and  extraordinary  modes  are  plotted  in  Fig.  6  as 
a  function  x,  where  x  —  —  r  cos  9  (negative  x  corresponds  to  the  weak  field  side).  From  the 
figure,  we  see  that  the  profiles  are  very  narrow  and  centered  around  x  -  0. 

A  typical  ray  trqectory  is  shown  in  Fig.  7,  where  we  plot  the  poloidal  plane  in  the  inner 
circle  and  the  toroidal  plane  (a  "top  view"  of  the  tokamak),  in  the  outer  circle. 
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Type  B  Tokamaks 

The  absorption  parameters  of  type  B  tokamaks  (e.g.  the  Princeton  Large  Torus)  are 
shown  in  Table  II.  The  dependence  on  9  and  9n  are  similar  to  that  for  type  A  tokamaks  and 
will  not  be  shown  here.  Here  we  are  interested  in  the  comparison  between  relativistic  and  non- 
relativistic  absorption  calculation,  and  the  effect  of  a  slight  divergence  of  the  incoming  waves, 
i.e.  (  ^  0  (see  Fig.  5  for  a  pictorial  definition  of  £). 

Since  all  published  calculations  on  electron  cyclotron  resonance  absorption  in  tokamak 
geometry  are  nonrelativistic,  it  is  of  interest  to  compare  nonrelativistic  and  relativistic  calcula¬ 
tions.  From  Table  II,  we  see  that  relativity  makes  only  a  small  difference  for  the  ordinary 
mode,  and  the  energy  is  essentially  100%  absorbed  near  the  center  in  both  cases.  However, 
there  is  a  big  difference  for  the  extraordinary  mode  propagating  at  ©  =  0°and  9„  =  180°.  The 
nonrelativistic  calculation  predicts  negligible  absorption  (because  of  the  shear  in  magnetic  field, 
kharhar  deviates  slightly  from  zero  as  the  mode  propagates  in,  and  the  absorption  is  roughly  1%). 
On  the  other  hand,  the  relativistic  calculation  for  the  same  parameters  shows  a  53%  absorption. 
This  is  due  mainly  to  the  broadening  of  the  cyclotron  resonance  through  the  energy  dependence 
of  the  electron  gyrofrequency.  For  0  -  30°  and  9„  -  125°,  we  see  that  the  extraordinary  wave 
is  fully  absorbed  in  both  the  relativistic  and  nonrelativistic  cases,  but  that  the  energy  deposition 
profiles  are  quite  different,  cf.  Fig.  8.  This  difference  is  enhanced  in  the  reactor  case  (type  C), 
as  we  shall  see  later. 

In  practice  the  heating  waves  will  typically  be  launched  from  a  horn  antenna.  Thus  there 
will  be  some  spectrum  of  launched  wavenumbers  (characterized  by  the  antenna  pattern).  The 
range  in  angle  of  deviation  |  of  the  launched  waves  could  be  as  large  as  ±  10°.  Table  II  shows 
the  results  for  f  —  0.5°,  5°  and  10°.  From  Fig.  9  and  Table  II  we  see  that  as  the  angular  diver- 
jence  increases,  the  mode  deposits  its  energy  farther  and  farther  away  from  the  center.  Thus, 
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as  a  result  of  the  spread  in  f,  the  composite  energy  deposition  profile  is  considerably 
broadened. 

The  transmission  profile  of  ordinary  and  extraordinary  modes  is  shown  in  Fig.  10.  Again 
we  see  that  the  profile  is  quite  narrow,  typically  about  5%  of  the  minor  radius.  Because  of  the 
high  temperature  (~  3  KeV),  the  extraordinary  mode  usually  does  not  deposit  its  energy  very 
close  to  the  center.  A  possible  remedy  for  this  is  to  tune  the  frequency  slightly  lower  than  the 
central  cyclotron  frequency.  Such  a  case  (o>  =  0.9  a)°(P)  is  shown  in  Figs.  10  (solid  curve)  and 
11.  We  see  that  the  maximum  energy  deposition  is  now  at  rmM! a  =  0.09  instead  of  0.33 
(Table  II). 

Type  C  Tokamaks 

The  absorption  characteristics  of  type  C  (reactor-size)  tokamaks  are  shown  in  Table  III. 
The  difference  between  relativistic  and  nonrelativistic  calculation  is  again  very  small  for  the  or¬ 
dinary  mode.  A  blow-up  of  the  transmission  curve  versus  x  is  given  in  Fig.  12.  For  the  ex¬ 
traordinary  mode,  however,  the  difference  between  relativistic  and  nonrelativistic  calculations  is 
significant.  The  energy  deposition  profiles  and  the  transmission  coefficient  profiles  for  the  ex¬ 
traordinary  mode  are  given  in  Fig.  13  and  Fig.  14.  The  effect  of  angular  divergence,  £,  is  rather 
strong.  We  show  in  Fig.  15  and  Fig.  16  the  results  for  (  —  0.5°,  5°  and  10°.  Again  it  is  evi¬ 
dent  that  as  (  increases,  (rmaJa)  also  becomes  larger. 

Finally,  we  study  what  happens  to  the  absorption  characterstics  if  we  use  electron  cyclo¬ 
tron  resonance  heating  to  heat  a  reactor  from  an  initial  temperature  of  2  KeV  (such  as  would 
result  from  ohmic  heating)  to  the  ignition  temperature,  say  7  KeV.  Blow-ups  of  the  energy 
deposition  curves  and  transmission  coefficient  profiles  are  shown  in  Fig.  17  and  Fig.  18.  We 
have  partly  optimized  the  heating  by  using  <o  —  0.95  t»'r.  It  is  clear  from  Fig.  17  that  the  ener- 
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gy  deposition  profile  becomes  broader  as  the  temperature  increases,  and  the  energy  is  deposited 
somewhat  further  from  the  center  of  the  torus  (rmax/a  =  0.02  at  Te  (0)  =  2  KeV  and  rmiJa  = 
0.28  at  7V(0)  =  7  KeV).  Note  that  the  launched  waves  are  the  same  for  Figs.  17  and  18,  and 
that  the  results  for  both  cases  indicate  wave  absorption  in  the  interior  of  the  device.  Thus 
these  results  indicate  the  possibility  of  heating  to  ignition  without  tuning  the  launched  wave  as 
the  temperature  rises. 

Discussion 

In  concluding  this  section  we  note  the  following  points: 

(1)  The  one  pass  absorption  is  usually  <  50%  for  small  tokamaks  (type  A).  However,  for 
larger  tokamaks  (type  B  and  type  C),  the  absorption  is  usually  essentially  100%. 

(2)  The  energy  absorption  is  usually  in  the  interior  of  the  device  and  can  be  further  centered 
by  tuning  the  frequency  slightly  lower  than  the  central  cyclotron  frequency. 

(3)  Relativistic  effects  on  the  absorption  coefficient  can  lead  to  significant  differences  (see  also 
Fidone  et  al.3). 

(4)  Angular  divergence  of  the  launched  waves  shift  the  energy  deposition  profile. 

(5)  Heating  of  a  tokamak  reactor  from  an  ohmically  heated  state  to  ignition  without  wave  tun¬ 
ing  is  feasible. 

(6)  Since  both  the  ordinary  mode  and  the  extraordinary  mode  are  fully  absorbed  for  type  B 
and  type  C  tokamaks,  the  choice  between  the  two  modes  should  depend  on  the  availabili¬ 
ty  of  space  (for  the  extraordinary  mode),  and  the  density  requirements  in  conjunction 
with  the  accessibility  condition.  In  particular,  if  we  are  interested  in  heating  high  density 
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tokamaks  [2  >  (at  0/w^)J  >1],  then  the  extraordinary  mode  is  the  only  possible  choice 
for  the  fundamental  harmonic. 

IV.  SCATTERING  BY  DENSITY  FLUCTUATIONS 

In  this  section  we  consider  the  effect  of  low  frequency  density  fluctuations  (such  as  those 
that  would  be  caused  by  drift  waves)  on  the  propagation  of  the  externally  launched  heating 
wave.  Strong  scattering  of  the  heating  wave  would  substantially  alter  the  picture  of  wave  propa¬ 
gation  and  energy  deposition  which  we  have  presented  in  the  proceeding  section. 

To  simplify  the  analysis  we  shall  consider  wave  propagation  perpendicular  to  B0  -  B0i„  in 
a  homogeneous  medium.  Since  the  frequency  and  parallel  wave  number  of  drift  waves  are  low 
compared  to  the  frequency  and  (perpendicular)  wave  number  of  the  heating  wave,  we  neglect 
them.  Hence,  on  scattering,  the  wave  frequency  is  conserved,  and  the  propagation  direction 
remains  perpendicular  to  B0.  Scattering  of  the  direction  of  propagation  in  the  plane  perpendicu¬ 
lar  to  B0,  however,  does  occur. 

First,  we  consider  an  ordinary  mode  heating  wave,  in  which  case  the  propagation  is 
described  by  the  following  wave  equation 

jvj?  +  -  ^y[i  +  a/KO/NlJf,  -  0.  (12) 

where  <op  is  the  plasma  frequency  in  the  absence  of  fluctuations;  N  is  the  average  electron  den¬ 
sity,  and  8 n  is  the  fluctuation  in  the  electron  density  which  depends  on  the  position  coordinates 
perpendicular  to  B„,  ij..  Fourier  transforming  (12)  with  respect  to  ra  we  have 

(— cJkJ+<uJ-«(?)£|l  *  (opN~'  £8flk  £k_k- 

k' 

A  A 

Replacing  <■>  by  «(0)  +  /-^,  where  [a»(0,)J  -  and  expanding  for  we  obtain 

the  usual  mode  coupling  equation 
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/  dEk/dt  -  N~'  £lF8flk£k_k..  03) 

k‘ 

where  the  coupling  coefficient  W  is  given  by 

W=<o}/(2w)  (14a) 

which  is  independent  of  k  and  k'  for  the  case  under  consideration.  A  similar,  although  algebra¬ 
ically  more  complicated  calculation,  also  yields  (13)  for  the  extraordinary  wave.  For  the  case  of 
the  extraordinary  wave  W  will  in  general  depend  on  k  and  k'.  However,  for  the  case  of  interest 
here  (  <■>  -  in  the  center  of  the  machine),  w  -  wce  throughout  the  device.  Adopting  this 
approximation  the  expression  for  W  simplifies  (cf.  Appendix  B  )  and  we  obtain 

W  =  [1  4-  (o>„ko,,)2)_1«pA2ci>),  (14b) 

for  the  extraordinary  mode,  which  is  independent  of  k  and  k'.  Applying  the  random  phase  ap¬ 
proximation  to  the  mode  coupling  equation,  Eq.  (13),  and  following  Ref.  11,  we  have 

HF  /»2i r 

-  2nk\ ^IVJo  dpSl2ksin(p/2)UF(4>  +  p)-F(4>)]  (15) 

where  F(d»)d<f>  is  the  energy  density  of  waves  with  wavenumber  vectors  whose  angle  to  the  x- 

X2ir 

F(<t>)d<t>  S(K± )  denotes 

the  two  dimensional  wavenumber  spectrum  of  the  low  frequency  density  fluctuations  (which 
are  perpendicular  to  B„),  defined  so  that  the  mean  square  relative  density  fluctuation  level  is 

<(8n/N)2>  —  S(K1)2vK±dK1, 

and  we  have  assumed  S  to  be  isotropic  (i.e.  it  depends  only  on  ^  -  |kJ).  We  note  that  the  in¬ 
tegral  operator  on  the  right  hand  side  of  (15)  has  eigenfunction  cos  pd>  (p  -  0,1,2 . )  and 

eigenvalues  -vr\ 

v„  -  4  irk  |  H'lV  sin2(pp/2)  S(2A  sin(/3/2))  dp  (16) 

As  discussed  in  Ref.  11,  a  reasonable  definition  of  the  scattering  length  can  be  based  on  the 

damping  (due  to  scattering)  of  the  lowest  order  angular  dependent  eigenfunction  ( p  —  1): 

/,  -  v>,.  (17) 

To  evaluate  v\  we  assume  that  S(k±)  is  Gaussian  (this  is  motivated  by  experiment,  cf.  Ref.  12) 
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SUX)  -  <  J^rJ  >  -yyj  exp  ( - 
and  il  «  k1,  in  which  case  (16)  becomes 

v's2T  1  Vi~  <<8*/A0J>.  (18) 

2  k2\K 

For  drift  waves,  reasonable  estimates  of  (0  and  <(8n/A02>  which  are  consistent  with  some 
recent  experiments  on  the  PLT  tokamak  are13 

-  pf1  and  <(6n/N)2>  -  ({„LB)-J.  (19) 

where  p,  is  the  ion  gyro-radius  and  L„  is  the  density  gradient  scale  length.  (Note,  however, 

that  these  estimates  do  not  apply  in  all  cases.  In  fact  very  large  relative  fluctuation  levels  have 

been  observed  on  the  Alcator  tokamak12.)  From  (14),  (17)  and  (18) 

/,  -  2(ir)-,/2t(*  v,)2/|  1F{21  (Lnlp,)2  Pi.  (20) 

For  6>  —  a>„  we  have 

k\g  s  o>„  (2  -  (tt>p/wc,)2]  ll  +  («„/«,) 2]-1  (21a) 

for  the  extraordinary  wave  and 

k  -  wf,ll  -  (w^/w,.,)2!  (21b) 

for  the  ordinary  wave.  For  example,  for  the  ordinary  wave  with  u2/<o2e  —  1/2,  (14)  (20)  and 

(21)  yield  l,  —  iL2/p ,  which  is  typically  very  large.  For  the  extraordinary  wave  with 

—  1  we  again  obtain  /,  —  8 L2/p,.  The  angular  scattering  of  the  propagation  direction  as 

the  wave  propagates  from  the  edge  to  the  center  of  the  tokamak  can  be  estimated  on  the  basis 

that  it  is  a  diffusion  process", 

A*  ~  (a/1,) 1/2  (22) 

where  a  is  the  tokamak  minor  radius.  Taking  L„  ~  a /2  we  obtain  to  be  in  the  range  Is  to 

3°  for  PLT  and  reactor  sized  tokamaks.  Thus  we  conclude  that  scattering  should  not  be  a  prob¬ 
lem,  if  the  characteristics  of  the  density  flucutations  are  approximately  consistent  with  (19). 
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V.  PARAMETRIC  INSTABILITIES 

In  this  section  we  consider  the  possibility  that  the  externally  launched  electron  cyclotron 
resonant  heating  wave  might  produce  parametric  instabilities,14  whereby  its  energy  will  feed  into 
lower  frequency  decay  waves.  Our  main  result  is  that  such  parametric  instabilties  are  unlikely 
under  conditions  that  will  apply  in  electron  cyclotron  resonance  heating  experiments  in  large 
tokamaks.  This  situation  is  in  contrast  with  lower  hybrid  heating  for  which  parametric  instabili¬ 
ties  are  thought  to  be  a  definite  problem.  The  reason  for  the  greater  stability  in  the  case  of 
electron  cyclotron  resonance  heating  can  be  crudely  explained  as  follows.  The  strength  of 
parametric  instabilities  depends  on  A,  the  displacement  of  a  particle  due  to  the  action  of  the 
pump  wave.  For  a  lower  hybrid  (///)  pump,  A"'  ~  E"'/(B0  o>//() ,  while  for  a  pump  with  frequen¬ 
cy  near  A“’  —  £,‘7(fl„  wir),  where  E„  denotes  the  heating  wave  (pump)  electric  field  am¬ 
plitude.  Furthermore,  for  a  given  input  power  flux  density,  P,  the  corresponding  pump  electric 
fields  scale  as  (£,f')2  ~  £/v“  and  ( £,‘0 2  —P/\"  where  vK,  is  the  radially  inward  component  of 
the  group  velocity.  Thus  for  fixed  P,  (A"'/A‘0  ~  (« !»„/«„,)  (v;,7v"')l/2.  Since 

—  (mjm,,) 1/2  and  v"'  »  we  see  that  (A "'/A"')  is  large.  Hence  lower  hybrid  waves 
might  be  expected  to  be  much  more  susceptable  to  parametric  instability  processes.  The 
preceeding  arguement  is  too  crude  to  be  taken  as  a  proof,  but  it  does  make  our  subsequent  con¬ 
clusions  plausible. 

If  a  window  is  used  to  isolate  the  wave  generation  region  from  the  tokamak,  then  the  heat 
loading  which  the  window  can  tolerate  sets  an  upper  bound  on  the  wave  energy  flux  density. 
For  continuous  operation  (i.e.  times  >  100  msec)  the  maximum  power  flux  density  which  a 
window  can  withstand  is15  —  10  kW/cm2.  Furthermore  we  note  that  this  flux  density  should  be 
quite  adequate  for  heating  a  reactor  plasma.  For  example,  using  numbers  from  the  UWMAK 
III  conceptual  reactor  design9,  we  find  that  with  an  incident  power  flux  density  of  S  kW/cm2 
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only  about  0.1%  of  the  available  surface  area  need  be  utilized  for  heating  ports.  Thus  if  we  find 
power  flux  density  instability  thresholds  which  greatly  exceed  10  kW/cm2  we  shall  conclude  that 
parametric  instabilities  are  unlikely  in  typical  reactor  situations. 


We  have  considered  parametric  instabilities  for  cases  in  which  the  decay  waves  are  elec¬ 
tromagnetic  and  for  cases  in  which  the  decay  waves  are  electrostatic.  We  found  that  elec¬ 
tromagnetic  decays  have  much  larger  instability  power  flux  density  thresholds  than  do  decays  to 
electrostatic  waves.  The  basic  reason  for  this  is  that  electromagnetic  waves  have  large  group 
velocities  (—  c),  and  hence  rapidly  convect  out  of  the  unstable  region.  We  therefore  limit  the 
subsequent  discussion  to  the  slower  moving  electrostatic  decay  waves. 

A.  Propagation  of  Electrostatic  Decay  Waves  in  a  Torus 


In  the  frequency  range  around  ioce  there  are  two  relevant  types  of  electrostatic  waves: 
cold  plasma  waves  (i.e.  the  upper  and  lower  hybrid  wave  branches)  and  the  electron  Bernstein 
waves.  (In  this  section  (<!>,  k)  denote  the  frequency  and  wavenumber  of  a  lower  sideband  elec¬ 
trostatic  decay  wave.)  The  cold  plasma  waves  are  characterized  by  k*pf  «  1  (where 
p,  -  vc/o>„  is  the  electron  Larmor  radius  and  vf  is  the  electron  thermal  speed),  and  have 
dispersion  relation 


which  can  be  solved  for  k± 


(23) 


k±  ~ 


2 

*  j 

Ota.  _  , 

.  “V  <*v 

2  1 

1  O  J  -2 

or 

3 

3 

1 

3 

Li 

#ii. 


from  which  we  see  that  propagation  ( k*  >  0)  occurs  in  the  ranges  u„i,  >  w  >  max(&>,,.  w,r) 
(the  upper  hybrid  branch)  and  min(<ufl,  «„)  >  «  >  o»w,  (the  lower  hybrid  or  Gould- 
Trivelpiece  branch)  where  a»2„  “  w2  +  «2,,  w2,  -  «£ll  +  (w(,/a»r(.)2l-\  and  we  have  made  use 
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of  the  approximation  me  «  m,.  We  take  the  heating  wave  frequency  to  match  at,.,.  in  the 
center  of  the  device.  Thus  if  the  magnitude  of  the  equilibrium  B-field  scales  approximately  as 
/?_l,  the  surface  where  the  decay  wave  frequency  matches  at„  must  be  at  R  >  Rn  (note  that 
the  decay  wave  frequency  must  be  less  than  the  pump  wave  frequency).  Thus  we  find  that  re¬ 
gions  of  propagation  of  cold  plasma  electrostatic  decay  waves  are  as  shown  in  Fig.  19  for  a  low 
density  case  (i.e.  at}  <  a>?e  everywhere),  and  in  Fig.  20  for  a  high  density  case  [a**,,  >  (w(")J, 
where  at^  and  w®,  are  atp  and  a>cc  evaluated  at  r  —  0],  Note  that  the  curves  at  -  represent 
reflection  surfaces  for  upper  hybrid  decay  waves;  the  curves  at  -  atp  represent  reflection  sur¬ 
faces  for  lower  hybrid  decay  waves;  and  the  curves  at  —  a>u/,  represent  cold  plasma  resonances 
of  the  upper  hybrid  decay  waves.  To  determine  what  happens  at  at  —  a>,.e  we  plot  the  dispersion 
relation  of  the  electron  Bernstein  mode  ( k{{  -  0)  in  Fig.  21  for  <u„„  <  2at,r  (i.e.  a <  3). 
We  note  that  the  Bernstein  wave  has  a  zero  in  kx  at  the  same  place  where  the  cold  plasma  wave 
has  a  resonance  (namely  at  -  atuh).  Since  Bernstein  waves  are  short  wavelength  modes 
(*iPf  ~  1  where  the  wave  is  in  the  center  of  the  region  >  at  >  at,.,.),  and  cold  plasma 

upper  hybrid  waves  are  long  wave  length  modes  (kx  p,.  «  1  in  the  center  of  the  propagation 
region),  we  see  that  with  finite  kH  and  Tc  the  cold  plasma  resonance  becomes  a  mode  conver¬ 
sion  surface  (cf.  Stix10).  Thus  an  upper  hybrid  mode  traveling  toward  the  at  -  at,,,,  surface 
mode  converts  to  an  electron  Bernstein  wave  which  then  travels  back  into  the  region  of  increas¬ 
ing  magnetic  field,  eventually  encountering  a  resonance  at  at  -  at,,,  (this  is  evident  from  Fig.  21 
which  shows  k  -*  »  as  at  -1 ►  «„).  With  regard  to  Figs.  19(c),  20(b)  and  20(c),  we  note  that 
the  Gould-Trivelpiece  waves  are  completely  enclosed  within  reflection  surfaces.  Thus  any 
Gould-Trivelpiece  decay  wave  is  forever  confined  within  the  high  density  region  of  the  torus. 
Such  decay  waves  could  experience  multiple  interactions  with  the  pump  as  they  bounce  around 
in  their  confined  region  and  transit  around  in  the  toroidal  direction. 
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B.  Parametric  Instability  and  Dispersion  Relation 


For  the  following  discussion  we  use  (Cl,  k)  to  denote  the  frequency  and  wavevector  of 
the  heating  wave  (pump)  and  ( a> ,  k)  for  the  frequency  and  wavevector  of  a  parametrically  dis- 
tabilized  low  frequency  perturbation,  «  «  ft.  The  parametric  dispersion  relation  for  electros¬ 
tatic  waves  in  a  magnetized  plasma  is  (cf.  Appendix  C  for  a  derivation) 


t(o),  k)  +  |k  •  A|2  «r(a»,k)[l+€,(w,k)l  {—  +  — |, 

i «+  * 


el 


(24) 


where  e±  =  e(o>±ft,  k±  ic),  «  ”  1  +  e,  +  th  and  «,  are  the  electron  and  ion  dielectric 
constants,  and  e±  is  defined  in  Appendix  C  (e±  —  1  in  the  dipole  approximation).  For  a  pump 
electric  field  E  -  E0  exp  (-iClt  +  iH  x)  +  (complex  conjugate),  the  quantity  A  given  in  Appen¬ 
dix  C  may  be  written 

A  -  -^-|z0(£K/ft)  +  *+lEJ(Cl+<*c)]  +  a_[£_/(ft-«jJ,  (25) 

where  E±  -  (2)-,/J  (£„?/  Eoy),  t±  -  (2)_l/2  (x0T/ y<>)  and  we  have  taken,  B0  -  B^0.  Note 
that  for  the  stated  definition  of  E0  the  temporally  averaged  pump  electric  field,  <L£>  — 
2E0  E0*.  As  a  point  of  reference  we  note  that  for  an  ordinary  wave  propagating  nearly  perpen¬ 
dicular  to  Bq  —  Bq10  k-k  x0  and  Cl  ~  <■>„,  we  have 


Eo“£o  *o> 

A  -  -c£§  z o/(B0Cl),  (26a) 

while  for  an  extraordinary  wave  with  k-k  x0  and  Cl  we  have 

E0-£o  (/* o  +  yo)/2'/2,  and 

A—c£ol— »  ft2/«2Xo  +  (ft2/w2— l)yol/(2,/2fl0ft)  .  (26b) 

Relating  Eg  and  £j  to  the  time  averaged  poynting  flux  (Pi,  we  have  from  Maxwell’s  VxE 
equation. 
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£o  -  (2 nil P/kc2)'/2  (ordinary  mode  )  ,  (27a) 

£o  “  ( 4it(IP/kc2)'/ 2  (extraordinary  mode)  .  (27b) 

The  two  relevant  types  of  parametric  decay  instabilities  which  we  consider  here  are  (1)  the 
nonresonant  decay  of  a  pump  (fl.ic)  to  a  lower  frequency  wave  (w-ft.k-*)  by  nonlinear  Lan¬ 
dau  damping  on  particles  satisfying  <w“»/r||V||;  and  (2)  three  wave  resonant  decay  of  the  pump 
(ft,)?)  to  a  low  frequency  wave  («,k),  with  n»<o,  and  a  lower  sideband  (w-fl.k-K). 
(Nonlinear  cyclotron  damping  at  a>  —  ^||V||  +  n<itie  is  not  considered). 

We  shall  first,  examine  nonresonant  decay  instability.  For  nonresonant  decay  instabilities 
the  wave  propagation  characteristics  of  the  decay  wave  determine  the  threshold.  For  decay 
waves  with  growth  rate  y  and  group  velocity  vK  we  may  define  k=y/vK,  the  inverse  exponentia¬ 
tion  length.  Then  the  total  number  of  exponentiations  of  the  decay  wave  as  it  convects  through 
the  unstable  region  occupied  by  the  pump  is  J k,dl ,  where  dl  is  a  differential  length  along  the 
ray  path  of  the  decay  wave.  A  criterion  for  convective  instability  is  then  that  the  total  number 
of  exponentiations  in  the  unstable  region  should  be  large.  On  the  other  hand,  for  decay  waves 
which  are  trapped  within  a  confined  region  (e.g.  the  region  within  tup>|w-n|  in  Fig.  20),  after 
a  sufficient  time,  the  decay  wave  sees  an  average  of  the  growth  rate  over  the  confined  volume. 
If  this  averaged  y  is  positive  (i.e.  the  averaged  parametric  growth  exceeds  collisional  and 
resonant  particle  damping),  then  an  absolute  instability  results.  Convective  and  absolute  non¬ 
resonant  instabilities  will  be  discussed  in  Secs.  VC  and  VD,  respectively. 

C.  Nonresonant  Decay  to  Electron  Bernstein  Waves 

For  nonresonant  convective  parametric  decay  the  most  dangerous  mode  is  the  electron 
Bernstein  mode  since  its  group  velocity  is  much  smaller  than  that  for  cold  plasma  waves.  This 
in  turn  results  in  a  larger  k,  for  Bernstein  waves.  Note  that  the  Bernstein  decay  wave  frequency 
is  less  than  the  frequency  of  the  pump  and  hence  only  propagates  on  the  low  magnetic  field 
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side  of  the  surface  where  the  pump  frequency  is  equal  to  the  cyclotron  frequency.  Since  the 
extraordinary  mode  must  be  launched  from  the  high  magnetic  field  side  of  the  torus,  it  will  al¬ 
ready  have  passed  through  its  region  of  strong  linear  absorption  before  it  can  excite  a  Bernstein 
wave  parametric  instability.  Thus  this  instability  is  primarily  of  interest  for  the  ordinary  wave. 


To  obtain  an  upper  limit  on  we  ignore  the  possiblity  that  the  Bernstein  wave  may 

convect  out  of  the  pump  region  as  it  propagates.  Thus  we  take  the  path  integral  from  the  upper 
hybrid  resonant  surface  £>  —  (more  accurately  the  mode  conversion  layer)  to  the  cyclotron 
resonance  layer  o»-a»rr  (cf.  Fig.  20),  where  «=0-Re(a>).  Since  k,  becomes  small  at  these  end 
points  we  may  roughly  estimate  f  k,dl  as  k°l0,  where  k°  is  k,  in  the  middle  of  the  region  where 
the  Bernstein  wave  propagates,  and  l0  is  the  length  of  the  ray  path  (roughly  the  distance 
between  the  &>—<■>„/,  surface  and  the  »-<■*„  surface).  An  upper  bound  on  /0  would  be  a,  the 
minor  radius.  To  estimate  X„  we  expand  the  parametric  dispersion  relation  for  small  k,  in  the 
nonlinear  Landau  damping  limit  and  obtain 


k,  -  - 


3«_ 


3*i 

+ 


-i 


Ml: 

klJ 


U+«,l2 


|l+e,  +  e,|2 


I  !+«/+«<■  1 2 


lm  (e,) 


(28) 


where  we  have  utilized  *»k  appropriate  for  Bernstein  waves  and  neglected  linear  damping 
processes.  For  nonlinear  Landau  damping  on  electrons  the  expression  for  k,  can  be  written 

-  (isrj  <»> 

where 


X-*iP,  , 

(  -  -|1+€,|2  |l+€/+€,|'2  Im(«,)  , 
Im(e,)-w  /0(X2)exp(-Xl)X":w2.ci),;2v<F0(v(,) 
Fo(v)“(2«,v2)"l/2exp-(v2/2v,2)  , 
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and 

v*«(n-w)//cl|  . 

Maximizing  f  we  have  £  —  0.35.  For  a  point  in  the  middle  of  the  unstable  region  A~1  and 
(d«_/3\)~l  for  Thus,  using  J  k;dlS  10  and  /„  ~  a,  we  obtain  the  following  esti¬ 

mate  of  the  instability  threshold 

|k-A|2  >30 (p„/a)  .  (30) 

Using  ft /k— c,  Cl~<occ—(oiv,  A~~(B0<oce)~t(2irP),/2,  and  Eqs.  (26b)  and  (27b),  with  X  —  1 ,  Eq. 
(30)  can  be  rewritten  as 


Po 


2xl04 


kW 


cm“ 


1  T' 

3/2 

_-l 

jlOkeF 

a 

(31) 


where  a  is  in  meters,  and  P„  is  the  threshold  power  flux  density.  (Recalling  that  our  estimates 
have  generally  favored  instability,  we  expect  P0  to  actually  be  underestimated  by  (31).)  Thus, 
if  as  previously  stated.  Pis  limited  to  about  10kW/cmJ,  we  see  that  P0  is  far  too  high  for  this 
parametric  process  to  occur. 


D.  Absolute  Parametric  Instability 


For  waves  which  are  trapped  within  the  interior  of  the  torus  [e  g.  Fig.  19(c)]  an  absolute 
instability  is  possible  if  the  average  parametric  growth  along  the  ray  path  exceeds  the  average 
collisional  and  resonant  particle  damping.  We  can  define  the  temporal  averaged  growth  rate  as 


<y>  =  lim  \t 


,-i 


J*o  y[r(r),  k(f)]  dt  , 


where  the  evolution  of  r  and  k  is  obtained  from  the  ray  equations.  Instability  occurs  if  <y>  is 
positive  (cf.  Appendix  D).  (We  shall  use  the  dipole  approximation  so  that  no  distinction 
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between  k  and  k  is  necessary.)  It  is  reasonable  to  assume  that  due  to  the  poloidal  dependence 
of  the  equilibrium  quantities,  there  is  no  constant  of  the  ray  motion  analogous  to  rk9  (which  is 
a  constant  of  the  motion  for  the  straight  cylindrical  case).  In  such  a  case  we  may  make  use  of 
the  Hamiltonian  form  of  ray  equations  to  infer  that  the  ray  wanders  ergodicly  over  all  regions 
of  phase  space  (i.e.  r-k  space  )  that  are  allowed  by  the  constants  of  the  motion  at  and  n  -  Rk+ 
(at  is  constant  because  the  equilibrium  is  time  independent  and  n  is  constant  because  of  toroidal 
symmetry).  In  fact  the  ergodic  wandering  of  rays  due  to  toricity  of  the  equilibrium  has  been 
explicitly  demonstrated  for  a  case  involving  lower  hybrid  waves  (cf.  Ref.  6  which  also  includes 
further  discussion  of  ray  ergodicity).  Thus  the  ergodic  theorem  may  be  utilized  to  convert 
<y>  from  a  time  integral  to  an  integral  over  phase  space  (i.e.  r,  k  space), 


J  y(r,k) 8[o»-a»(r,k)] 6lk+-(n/R)]drdk 

<y>  -  J — j. - ,  (32) 

J  Sldt-dt(r.k)]Slk^-(n/R)]drdk 

where  at—  fl-Re(w)  is  the  decay  wave  frequency  with  w(r,  k)  a  solution  of  (23),  and  8  (. . .) 


denotes  the  delta  function. 


Again  making  the  electron  nonlinear  Landau  damping  and  the  dipole  approximation  in 
(24)  we  obtain  the  temporal  growth  rate 

[|k-A|Jf-/m(c_)I  .  (33) 

The  term  |k-A|Jf  represents  parametric  growth  [cf.  Eq.  (29)  for  the  definition  of  f)  and  is 
counteracted  by  the  linear  damping  term  /m(c_).  For  /m(e_)  we  utilize  the  collisonal  contribu¬ 
tion  only,  since  the  resonant  particle  damping  at  a»— /cMvn  may  be  made  small  by  making 
at/kn  »  v,: 


where  v  is  the  electron-ion  collision  frequency,  and  v/at  «  1  is  assumed.  From  (33)  and  the 


i 
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assumed  ergodicity  of  the  decay  wave  ray  trajectory,  we  obtain  the  instability  condition 


<|k-A|2f>  5  <Im(«_)>  ,  (34) 

which  follows  from  <y>  >0.  Again  we  utilize  rough  estimates  to  obtain  an  idea  of  whether 

the  existence  of  such  an  instability  is  plausible.  Taking  0.3S  (its  maximum  value), 

A2~(B0<oce)~I2irPc>  (o?e~<02e  and  Im(e_)  ~  we  obtain  the  instability  condition 

<P>  >  5J(.N  T,  c)(vlu>f)(kp,)~2  ,  which  with,  v  —  3  x  10-6  N X7V*3/2  (where  N  is  in  cm-3, 

T,  is  in  eV,  and  X  is  the  usual  Coulomb  logarithm  A — 15  )  yields  the  threshold  power  flux  aver* 

aged  over  the  plasma  volume 


<P>  o 


2.2 


kW 

cm2 


10  keV  1/2 

N 

i/2 

0.1 

Te 

10ucm-3 

Recalling  that  k2p 2  «  1  for  the  validity  of  (23)  we  see  that  k2p}  —  0.1  is  a  reasonable  upper 
bound.  Also  assuming  that  heating  ports  take  no  more  than  one  percent  of  the  tokamak  wall 
area,  we  have  P  >  10 2<P>.  Thus  P0,  the  threshold  power  flux  density,  satisfies 


P0  >  (220  kW/cm2)  (lOkeV/T^lN/dO^cm-3)]2'1 ,  (35) 

and  again  we  see  that  P0  is  too  large  for  instability  to  occur  in  typical  situations  for  a  tokamak 

reactor.  In  extreme  conditions  for  a  present-day  lower  density  device  the  instability  might  be 

possible.  For  example,  with  P—5  kW/cm J,  T,  —  IkeV,  instability  occurs  for  N  ^  4  x  101J  cm-3. 


E.  Resonant  Decay  to  an  Ion  Cyclotron 
Wave  and  a  Cold  Plasma  Wave 


The  wave  number  of  the  pump  satisfies  *c  —  The  frequency  of  the  cold  plasma  de¬ 
cay  wave  is  slightly  (by  about  <uc/)  below  the  pump  frequency;  so  it  too  is  of  order  Since 
the  cold  plasma  wave  is  electrostatic,  we  assume  that  its  phase  velocity  is  much  less  than  c,  i.e. 

»«/*  «  c  ~ 

Thus  we  have  k  »  k.  This  justifies  the  use  of  dipole  approximation.  Furthermore,  we  can 


write 


kpi  ~  wceVi/(ucic)  -  1.9 (»»///»,) 1/2  T,U2(keV)  , 
where  mp  is  the  mass  of  a  proton.  Since  k  »  k,  we  have 

Jtp,  »  1.9  (mJmpYn  T^HkeV)  . 

Thus,  for  T,  >  200  m,/mp  cK,  we  have 

k2p?  »  1  , 

and  e~l+<(+e,  for  the  ion  cyclotron  wave  can  be  written 

r-  1 +k~2\jg  -  V2  uci/[(2ir)U2kipi  v2(o>— «f/)]  .  (36) 

Since  the  decay  occurs  for  ion  cyclotron  waves  which  have  k2  «  k2  and,  since  the  di¬ 
pole  approximation  is  valid,  the  wavenumber  of  the  high  frequency  decay  wave  is  also  almost 
perpendicular  to  B.  This  implies  that  the  high  frequency  decay  wave  is  an  upper  hybrid  wave 
existing  near  the  uppelr  hybrid  resonance  surface,  which  is  toward  the  plasma  edge  on  the  low 
magnetic  field  side  of  the  tokamak.  Thus,  this  decay  process  is  unimportant  for  an  extraordi¬ 
nary  wave  which  must  be  launched  from  the  high  magnetic  field  side  and  propagate  through  the 
resonant  absorption  region  (o»  =  &>„)  before  it  reaches  the  region  where  this  process  can  occur. 
Hence,  we  only  consider  this  process  for  an  ordinary  wave. 

For  a  three  wave  decay  instability,  the  parametric  dispersion  relation  [Eq.  (24)1  in  the  di¬ 
pole  approximation  reduces  to 

*«_  +  |k-A|2«,(l+C/)  -  0  ,  (37) 

where  <_  is  given  by  (for  k2p2  «  1) 

“  1+  +  a H  Ik2,  (38) 

with  ej.  -  -w£/(«2-«£)  and  cn  -  -u^/u2.  By  Taylor  expansion  of  Eq.  (37),  we  obtain 


where  y  —  Im(«). 


££r>_.  miv.  (l«0. 


(39) 


nmm 
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The  quantity  of  interest  is 


=  y2/|  '/gx'/fx-  -4-  I*x  U)-  Acx_  (x)]|  , 


=  r  /  I  ’xx'pr-  ^ 


(40) 


where  the  minus  subscript  refers  to  the  lower  sideband  (upper  hybrid)  wave.  In  terms  of  y.  the 
instability  threshold  criterion  is  ij  >  1  (cf.  Rosenbluth16).  Inserting  the  expression  for  y2  ob¬ 
tained  from  Eq.  (39)  into  Eq.  (40),  and  applying  the  relations  v,*  -  Mbkx 
— 0«/8*j.)  Oc/Ow)"1  kx/k±,  «,  -  -1-c,  -k2ks},  we  obtain 

I*!/*,!  Ml2 


*4X* 


From  Eqs.  (36)  and  (38),  we  obtain 


4  I 

8*1  8*i  dx 


(41) 


8«_/8*i  -  -  (l+«„)/*4  . 

df/8^  -  (Ac^Xi,)-'. 


rf«i/dx 


77  1 


(42) 

(43) 

(44) 


where 


dtjdx  —  —u)~2(dcj^/dx  +  dotlt/dx)  . 


(45) 


We  observe  from  Eqs.  (43)  and  (44)  that  d(kX-~k^ldx cannot  be  made  zero  in  the  region  of 
interest  (i.e.  the  low  magnetic  field  side  of  the  tokamak),  where  all  the  terms  on  the  right  hand 
side  of  Eq.  (44)  have  the  same  sign.  We  also  observe  from  Eqs.  (41)-(44)  that  y  maximizes  at 
fip»n  *„  so  that  the  second  term  on  the  right  hand  side  of  Eq.  (44)  can  be  neglected.  Thus, 
combining  Eqs.  (41)-(44),  we  obtain  the  following  threshold  criterion  for  the  instability 

To  evaluate  y  numerically,  we  note  that,  for  small  angles  (•,)  between  k  and  the  magnetic 
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field,  |k-A|2  -  fc2A2  sin2#,  -  fc2A2#2  for  an  ordinary  wave.  From  Eq.  (45),  we  can  approxi¬ 
mate  dtjdxby  a-1,  where  a  is  the  plasma  radius.  Thus,  Eq.  (46)  can  be  expressed 

V  -  |-^|  sin2#,  £  ka  >  1  (47) 

K  ««  p; 

Since  the  wave  frequency  is  approximately  the  upper  hybrid  frequency  at  the  instability  point 
and  the  cyclotron  frequency  in  the  center,  we  have  +  «2,(1— r/R)*'  «£.  Hence 
2r//?  and  Eq.  (47)  can  be  written 

V  “  2|-^-|  # ?-jr  Af  ka  >1  .  (48) 

"  Pe 

Substituting  A2  -  2ir  P  c/B2wit  into  Eq.  (48),  we  obtain  the  instability  condition. 


,  -  0.2*|  — I— fl(„  T„(key)  >  1  ■ 

For  example,  with  kp ,  <  1,  r/R  <  0.2,  #/  <  0.3  (about  10®),  Eq.  (49)  yields 


P  >  (200  kW/cm2)  \kx/k±\  B0  Te3,2a~l.  (50) 

with  B0,  T,  and  a  in  tesla,  keV  and  meters,  respectively.  Equation  (50)  shows  that,  except  for 

\kjkx\  ?  10-2  (side-scattering),  this  instability  would  not  take  place.  The  case  of  side¬ 
scattering  requires  a  more  involved  analysis  along  the  lines  of  Ref.  17  which  we  defer  for  furth¬ 
er  study. 


VI.  CONCLUSIONS 


Neither  parametric  instability  nor  scattering  by  existing  density  fluctuations  appear  to  be 
serious  problems  for  heating  tokamaks  by  the  electron  cyclotron  resonance.  The  accessibility 
conditions  do  impose  serious  limitations  on  the  maximum  density  that  can  be  heated.  Howev¬ 
er,  beta  limitations  on  stability  may  result  in  tokamak  parameters  consistent  with  the  accessibili¬ 
ty  requirements.  (The  acessibility  condition  is  further  relaxed  for  2  heating,  which  we  shall 
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discuss  in  a  future  publications21.)  The  predictions  of  a  linear  ray  trajectory  code  have  shown 
that  for  large  tokamaks  the  wave  energy  absorption  is  essentially  complete  on  one  pass,  and  that 
this  absorption  is  in  the  plasma  interior.  Furthermore,  it  appears  feasible  to  heat  a  reactor  plas¬ 
ma  from  an  ohmically  heated  state  to  ignition  without  wave  tuning.  In  addition,  by  properly 
launching  the  heating  wave,  it  may  be  possible  to  attain  a  degree  of  control  over  the  energy 
deposition  profile  (this  may  be  desirable  in  order  to  eliminate  desruptions  due  to  internal 
kinks).  We  emphasize  that  these  linear  predictions  should  be  accurate  in  the  light  of  our  tenta¬ 
tive  findings  that  parametric  instabilities  and  scattering  are  unimportant.  Thus  we  conclude 
that,  subject  to  uncertainties  in  the  lutimate  development  of  gyrotron  wave  sources  and  in  the 
stable  beta  limit,  electron  cyclotron  resonance  heating  has  many  favorable  features  which  make 
it  attractive  for  tokamak  applications. 
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Appendix  A 

THE  LINEAR  DISPERSION  RELATION 


The  conditon  for  electron  cyclotron  resonance  is 


It  can  be  shown  that  if 


10  -  *„v„  -  <o„/y0  -  0  . 


(A.l) 


k»c/w„  <  \Jc  ,  (A.2) 

where  vr  is  the  electron  thermal  speed,  the  third  term  in  (A.l)  is  a  more  sensitive  function  of 
electron  thermal  spread  than  the  second  term.  [Equation  (A.2)  can  also  be  written, 
6  <  2.7  T#1/J'  with  T,  in  keV  where  sin  0  —  knc/<o.\  Since  variation  of  the  third  term  is 
caused  by  the  relativistic  factor  y,  it  is  expected  that  under  condition  (A.2)  the  relativisitc  effect 
will  play  an  important  role  in  the  wave-electron  resonance.  Therefore,  to  obtain  a  general  for¬ 
mulation,  one  needs  to  employ  the  relativistic  Vlasov  equation  and  the  full  set  of  Maxwell 
equations.  Letting  k  —  kx  tx  +  ku  en,  and  assuming  a  Maxwellian  electron  distribution, 

fo  -  (2ir)~i,29;3  exp[-p2/202]  . 

where  p  is  the  momentum  and  0,  =  m,v,,  we  have  derived  the  following  dispersion  relation 

lO 

after  some  standard  algebraic  manipulations  . 

D  -  det|</„|  -  0,  (A.3) 

where 

d„  —  1  -fcjfcVa*2  +  <  J}M>/9}k±  w  , 

H 

d„-\ ~  kHeW  +  £  m}  <  P±J?M  >  /$}-  . 
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dn  “  1  -  XjV/u2  +  £  W  2  <  />?i^2(X)  >  /9  2<o  , 

n 

idxy  “  ~ id yx  “  ^  nmc<Op(o„  <  p^-/B(X)y,  (X)  >  /0 2  ^  o>  , 

« 

da  “  4*  “  k±kHc2/(o2  +  £  nmeo>2u„  <  p„72(X)  >  /» 2Xa  a*  , 

n 

idyz—idp  --'£*>}<  pj>\\  J„M  y;(x)  >  /»2«  , 

n 

X  -  kj>Jmto)„  . 

/;(x)  -  dj„M/dx , 

<a>  s  2.  /; pA £  „  ■ 


(A.4) 


To  evaluate  integrals  of  the  form  (A.4),  we  assume  X«1  and  expand  the  Bessel  funtions  in 
terms  of  X.  Note  that  Landau  path  is  implicit  in  the  Pn-integration.  To  lowest  order  in  X,  we 
obtain  the  cold  plasma  dispersion  relation  [Eq.  (1)  of  the  main  test).  To  next  order  in  X,  there 
are  two  types  of  contributions  to  the  integral.  The  resonant  electrons  give  rise  to  the  imaginary 
part  of  <<?>  (through  the  residue  of  the  poles)  and  consequently  cause  the  wave  to  damp. 
The  nonresonant  electrons  only  add  an  insignificant  contribution  to  the  real  part  of  <G>.  For 
the  present  study,  we  will  neglect  these  (higher  order)  nonresonant  electron  contributions. 
After  some  algebra,  we  obtain  the  final  expressions  for  the  du' s, 

d„  “  1  -  *iiC2/«2  -  w2/(«2-tt£)+/  (2 ir)l/2«2m,c  /( O)/80e5<o  , 
dyy  -  1  -  k2c2/<o2—<Op/(<o2—ot2e)  +  /( 2ir)l/2  a>2m,c/(O)/80’a»  , 
da  -  1  -  k2c2/w2—Up/<o2  +  i(2nV,2<o}k2c  /(2)/80*«m,ci>2,  , 
idxy  -  ~idyx  -  -a»„ai2/w(w2-<u2,)  +  i(2ir)ulu2mec  /( O)/80,5a> 
da“  da-  *||*j.c2/«2+/(2ir)I/J  w2  /(l)/80/a>  w„  , 
id,,  -  -  -/(2ir),/2»**1c  /(1)/80,5«u«„  , 

where 

/(«)  “  J^"  PiFis.pJ  dp±  , 

F(s,pj)  -  ( ( 1  +px2  /m,2c2)  (X|?c J-«2)  +«£) -l/2 

x{y_p,?_exp[-(pJ.2+pl?_)/202l+dy+F|f4expl-(p12+plf+)/202]), 

^n±  “  («2-*ifc2)_l(*HC«„±«I(l+p12/m#2c2)(*|fc2-»2)+«i)), 

y±  -  (l+(Pi2 +Pi?±)/m,2c2]  , 
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Re  {*,cl«i/U»J-*i?c2)-l]1/J)  )  if 

00  ,  if  <o/kn<c. 


1  ,  if  t»/k„>c  , 
0  ,  if  »/*n<c  . 


Appendix  B 

CALCULATION  OF  W  FOB  THE  EXTRAORDINARY  MODE 
From  Maxwell’s  equations  we  have 

VxVx  E  -  Ke  - 


c1  - 


where  m  is  the  mobility  tensor 


M  “  cw„  [BoWe-w2)]-1  1/  («-»„)  *oxl  A 

where  /  is  the  unit  tensor  and  we  have  assumed  that  E  is  perpendicular  to  Bo  and  has  time 
dependence  exp(-iur).  Fourier  transforming  in  space  we  have 


-kxkx/-^  /  + 

-  C1  -  C1  — 


.  Ek  -  -  •  £  */ik_k.  Ek'.  (B-l) 

C2  -  V 


With  8/»— *0,  Eq.  (B-l)  yields  the  well-known  properties  of  the  linear  extraordinary  mode  pro¬ 
pagating  in  a  homogeneous  medium:  the  dispersion  relation  is 

D(k.m)  -  *2-4+4  -  0. 

c2  c2 

and  the  polarization  vector  (  a  unit  vector  in  the  direction  of  E)  is 

to  -  (l+p2)-,/i  [y0+/pxo]  , 
for  wave  propagation  in  the  x  direction,  and 
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Note  that  for  <u  -  p-  1  and  e0  -  (2)-1/2  (y0-Hx).  Let  ek  denote  the  polarization  vector  for 
a  wave  propagating  in  some  arbitrary  direction  k  -  (k  cos/8,  k  sin/9),  then  ek  —  R-eo  where  R 
is  the  rotation  matrix 


R 


cos/3  -sin/3 

sin/3  cos/3 


Dotting  (B-l)  with  ek*we  obtain 

D(k, at)  £k  -  -(l+p2)4 iriuec'1  £  ek*  p  ek  8nk_k£k  .  (B-2) 

k’  " 


where  we  have  used  Ek  =  ek  £k.  As  for  the  ordinary  mode  case,  we  replace  <■>  by  «(0)  +  ib/bt, 
2>(«<0\  k)  s  0,  and  expand  (B-2)  for  «<0)  »  d/dt.  The  result  is  the  mode  coupling  equa¬ 
tion,  Eq.  (13),  with  IF now  given  by 


W  -  —  O+P2)  Airio>c~2Ne  ekM  ek 

For  u  —  o>„  this  expression  for  W  yields  (after  some  algebra) 


-l 


(B-3) 


[1+^t]  [cos^  "^)  +  > 0  -/3)]  . 


(B-4) 


which  gives  Eq.  (14b). 


Appendix  C 

PARAMETRIC  DISPERSION  RELATION  FOR  ELECTROSTATIC 
DECAY  WAVES  IN  A  MAGNETIZED  PLASMA 

In  this  appendix  we  use  the  ponderomotive  potential  method19  together  with  the  expres- 
sion  for  the  ponderomotive  force  in  a  magnetized  plasma*u  to  derive  the  parametric  dispersion 
relation  for  electrostatic  waves.  (  The  dipole  approximation  is  not  necessary  and  is  not  used  in 
what  follows.) 

Let  <f>±  and  8/t±  denote  high  frequency  potential  and  second  order  parts  of  the  electron 
density  perturbations  at  (a»±ft  and  k  ±  *  ).  Let  d>'  and  bn‘t  be  low  frequency  potential  and 
electron  density  perturbation  at  («,  k  ).  Let  E0  denote  the  complex  pump  electric  field  ampli¬ 
tude  such  that  the  actual  pump  field  is  E<>  exp  (— iClt+in  •  x)  +  (complex  coi\jugate). 

a.  High  frequency  components:  From  Poisson’s  equation  and  the  definition  of  the  plasma 
dielectric  constant 

(k±*c)J  «±  d>*  --4irAfe  («»*/*)  (C-l) 

where  c±  -  <(o>±n,k±K>  are  the  dielectric  functions  for  waves  with  frequencies  m±Cl  and 
wavenumber  k±K.  From  the  continuity  equation  and  the  cold  plasma  electron  momentum 
equation  we  obtain 

^  X(O)  •  E. 

4n(v  ±  (Dein*  ±  I  ±Hd  •  x *(fl)  .  g  *  *"*  “  °*  (C-2) 
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where 

/ft«i)„(ft2-a>c,)_1  o 

0  1 

From  (C-l)  and  (C-2)  we  obtain 


nJ(ftJ-a>2?)-1 
x  (n)  —  /Ow„(ft2-<u«)-1 

m 

0 


x(n)E0 


(k±K  Vt***  -  /  (k+K)  x*(ft)  Eo* 


b.  Ponderomotive  potential  The  ponderomotive  force  [at  (w,  k)  ]  on  an  electron  in  a  mag¬ 
netic  field  in  the  presence  of  high  frequency  electric  fields  E0,  and  [-/(ki*)^*]  is  ik  <t>p  e 
where 


<t>.  -  — — — r  {E0*-x,['(k+*)^+)  +  I/(k-u)*-l-X*Eo) 
me£ll  -  ~ 

Inserting  (C-3)  into  this  expression  for  4>p  we  obtain 

4 wk2ed>p  -  -|k-A|2  e'finl  (C-4) 

where  e±  -  |[(k±K)A]  (JrA)-I|.  and  A  -  -ex  E 0/(meft2),  and  we  have  made  use  of  the  her- 

mitian  property  of  x 


c.  Low  frequency  perturbations.  The  electron  and  ion  dielectric  constants  at  (o>,k)  (denot¬ 
ed  a,  and  a,  )  give  the  density  response  to  forces  applied  on  the  electronic  and  ionic  plasma 
components,  thus, 

k2et  (^'+0,)  “  4 irebn1',  (C-5a) 

k2e,  d>'  ”  —4 irebnl.  (C-5b) 

(+p  acts  only  on  the  electrons  since  (m,/m,)«  I.)  In  addition  we  have  Poisson’s  equation  for 


36 


NRL  MEMORANDUM  REPORT  4028 


the  low  frequency  components 

*V  =  4w(>(8m/-8m/)  .  (C-6) 

Eliminating  <ftp  </»'  8m,!  and  8m/  from  (C-4)-(C-6)  we  obtain  a  general  parametric  desperison  re¬ 
lation  for  electrostatic  waves  in  an  applied  magnetic  field 

e  +  | k-A| 2  e,,(l+«,)  Up +  /«+)+(? -/«-)]  =  0.  (C-7) 

where 


€  =  I  +<,.  (o>,  k)  +  €,  (o).k)  . 


In  the  dipole  approximation  (k  «  Ar),  we  have  e±  -  1. 


Appendix  D 

VALIDITY  OF  THE  ERGODIC  INSTABILITY  CONDITION,  <y>  >  0 


The  exponentiation  of  wave  energy  is  given  by  exp[2j^ y(r  )«/r’l,  where  the  i  dependence 
of  y  is  due  to  the  lime  dependence  of  r  and  k,  which  evolve  according  to  the  ray  equations. 
Defining  y  =*  <y>+Ay  we  can  write 

<exp[2^  y(i')di']>  =  exp(2<y>/  +2  <(J^  Ay (/')  <//']2>),  (D-l) 

which  would  be  exactly  true  if  J  \y(t')dt'  were  a  Gaussian  random  variable  (which  it  is  not). 
In  order  for  <y>  >  0  to  be  a  reasonable  instability  criterion,  the  second  term  in  the  exponent 
of  (D-l)  should  be  small  compared  to  the  first  term,  or 

<y >1  »  <[J(Ay (i')dl']2>  .  (D-2) 

Defining  a  correlation  function  C’(.v)  =  <Ay(f+.v)Ay(f)  >,  and  a  correlation  time  t,,  such  that 

<(Ay)2>  r,  -2f~C(s)ds  . 

we  can  write  (D-2)  as 

<y>  »  <(Ay)J>  t,  ,  (D-3) 

for  t  »  t,  .  For  the  strongly  crgodic  case,  r,  is  of  the  order  of  a  radial  transite  lime  of  the 
ray.  Also,  assuming  <(Ay)2>  ~  ( < y > ) 2,  (D-3)  becomes 

I  »  <y>  r,  . 

Thus  the  average  growth  on  one  radial  transit  must  be  small  in  order  for  <y>  >  0  to  be  a 
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Table  I.  Absorption  characteristics  of 
Type  A  Tokamaks  ("small  size") 


MODE 


ABSORPTION 

(%) 

'mJo 

11 

0.03 

11 

0.03 

10 

0.03 

10 

0.03 

19 

0.03 

28 

0.02 

39 

0.02 

19 

0.16 

46 

0.02 

0.8 


5°  |  10°  0.8°  0.8 


0°  0.8° 


10°  0°  0.8°  I  0.8 


1.2 


30°  I  180°  0.8° 


30°  180°  0.8°  0.8 


30°  135°  0.8°  1.2 


40°  180°  0.8° 
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I'i(  4  —  Minimum  accessible  value  of  (w^/ co in  the  center  of  Ihc  plasma  versus  .V1  for  (a)  the 
extraordinary  mode  and  (h)  Ihc  ordinary  mode. 
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Fig.  6  —  Transmission  coefficient  for  Ihc  ordinary  and  extraordinary  modes  in  a  type  A  lokamak  versus  v.  The  long 

dash  curve  ( - )  is  for  Ihc  ordinary  mode  with  “  0  8,  «  -  5°,  H0  -  0°.  The  extraordinary  mode  is  given 

by  Ihc  short  dash  curve  (— )  for  u>ra “  1.2,  the  solid  curve  ( - )  for  u>r0/o>«.  “  1.0  and  the  dash  dot  curve 

(—  •  — )  for  u>p0/ui".  “  0.8.  For  Ihc  three  extraordinary  mode  eases  H  -  30°,  »n  —  180°  In  all  eases  B0  -  12  KG 
and  Tr  -  350  eV. 


Fig.  7  —  An  example  of  a  ray  trajectory  for  a  type  A  lokamak  for  the  cxlradinary  mode  with  M  -  30°.  Hn  -  135 
“  I-2*  Te(0)  —  350  cV  and  Bq  —  12  KG.  The  inner  trace  is  in  the  po loid.il  plane  and  the  outer  trace 
in  the  toroidal  plane. 


—  lincrgy  deposition  versus  radius  of  magnetic  surface  /  for  a  type  H  lokamak.  (a)  relativistic  and 
(hi  non  relativist iv.  H  -  JO*.  «„  -  125*.  —  0.81.  r. (01  -  2‘>  KeV  and  H„  -  52  K(i. 


ENERGY  DEPOSITION  (ARBITRARY  UNITS) 
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I'ig  V  —  Cncrgy  deposition  versus  radius  of  magnetic  surfaec  r  for  a  type  B  loKamaK,  (a)  (  -  0  5"  (h) 
f  -  5”  and  (el  f  -  10’  H  -  JO*.  «„  -  125".  -  O  KI.  /;.<())  -  2  4  KeV  and  H„  -  52  K(i 
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Fig  10  —  Transmission  coefficient  of  ordinary  and  extraordinary  modes  in  a  type  B  tokamak  versus  v  -  -rcosW. 

The  dash  curve  ( - — )  is  for  the  ordinary  mode  with  top  ™  0.81,  H  -  2°  and  -  0°.  The  extraordinary 

mode  is  the  dash  dot  curve  I—  '  —I  with  urn “  0.81,  M  -  .10”,  «  -  125°;  the  solid  curve  has  the  same  parame¬ 
ters  but  with  m  -  0.9  oi®..  H(,  -  32  KCi  and  7V(0)  -  2.89  KeV 
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Fig.  1 1  —  Energy  deposition  versus  radius  of  magnetic  surface  r,  (a)  u  -  tu®. .  <b)  u  -  O  f  «®..  a>,o/«,r  “ 

0  *1.  «  -  30*.  #n  -  125*.  H  -  32  KG,  r,.<0>  -  2  89  KeV 
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Fig.  17  —  F.ncrgy  deposition  versus  r  for  .1  type  C  lokamak,  (ill  r,.(0)  -  2  KeV  and  <h)  rr(OI  * 
7  KeV  H  -  10*.  (*0  -  135°,  Hn  -  40.5  KG  and  Wo.”  -  0.95. 
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Fig.  19  —  (a)  Upper  hybrid  frequency  electron  cyclotron  frequency  to»l(1),  and  plasma 

frequency  versus  position,  v  “  R0  -  R.  where  R  is  the  major  radius  and  R0  is  the  major 
radius  at  the  center  of  the  minor  cross-section  (r  —  0).  <u,  and  u>j  arc  frequencies  of  possible 
decay  waves  <W| .  wj  <  The  density  of  low,  i.e  u>p  <  oj„.  everywhere  (b)  Propagation 

region  for  a  cold  plasma  decay  wave  with  frequency  w,  (an  upper  hybrid  wave),  (fi  is  the 
pump  frequency.)  (c)  Propagation  region  for  a  cold  plasma  decay  wave  with  frequency  iij  (a 
lower  hybrid  wave). 
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l'8  10  (continued)  —  (a)  Upper  hybrid  frequency  (<»„„>,  electron  cyclotron  frequency  (i«,,.).  and  plasma  fre¬ 
quency  l«in)  versus  position.  v  —  Rn  -  N.  where  K  is  the  major  radius  and  K„  is  the  major  radius  at  the  ccrucr 
ol  the  minor  cross-section  (r  —  0).  <".» ,  and  are  frequencies  of  possible  decay  waves  ( i*j| ,  <,> j  <  The  den¬ 

sity  of  low,  i.c.  mp  <  miv  everywhere,  (h)  Propagation  region  Tor  a  cold  plasma  decay  wave  with  frequency  (u I 
(an  upper  hybrid  wave).  (It  is  the  pump  frequency.)  (c)  Propagation  region  for  a  cold  plasma  decay  wave  with 
frequency  ill j  (a  lower  hybrid  wave) 
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